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We investigate the von Neumann entropy of a block of subsystem for the valence-bond solid (VBS) 
state with general open boundary conditions. We show that the effect of the boundary on the von 
Neumann entropy decays exponentially fast in the distance between the subsystem considered and 
the boundary sites. Further, we show that as the size of the subsystem increases, its von Neumann 
entropy exponentially approaches the summation of the von Neumann entropies of the two ends, 
the exponent being related to the size. In contrast to critical systems, where boundary effects to 
the von Neumann entropy decay slowly, the boundary effects in a VBS, a non-critical system, decay 
very quickly. We also study the entanglement between two spins. Curiously, while the boundary 
operators decrease the von Neumann entropy of L spins, they increase the entanglement between 
two spins. 
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I. INTRODUCTION 

Recently much research has been undertaken to un- 
derstand the subtle interplay between quantum entangle- 
ment and quantum criticalitu for spin systems [H [Ij H S 
Si, 0, S,[l- Vidal et al. [|] showed that the entangle- 
ment between a block of contiguous spins and its comple- 
ment in the ground state of the Ising model shows differ- 
ent behaviours for the gapped and gapless cases (critical 
and non-critical). The entanglement of the VBS ground 
state of the much-studied Affleck-Kennedy-Lieb-Tasaki 
(AKLT) model 0, Hil, [ll is considered in Refs. SQ, 
and very recently Campos Venuti et al. studied this 
state's long-distance entanglement property The 
entanglement of the fermionic system was studied in 
Refs. [il, [13. 

While much theoretical work in this area has focused 
on periodic boundary conditions, the open boundary con- 
dition has also attracted recent attention 15, Laflo- 
rencie et al. [15] numerically studied the boundary effects 
in the critical scaling of entanglement entropy (von Neu- 
mann entropy of a block of spins) for the (gapless) ID 
XXZ model, and found the entanglement entropy slowly 
decays away from the boundary with a power-law. This 
result can be interpreted as stating that in critical sys- 
tems, the boundary effects to the entanglement entropy 
is quasi-long-ranged; i.e. there is a quasi- long-ranged en- 
tanglement between the boundaries and the subsystem 
in question. This agrees with the fact that entangle- 
ment entropy increases logarithmically with the size of 
the subsystem in critical systems By contrast, 

the entanglement entropy for non-critical systems satu- 
rates to a constant bound when the subsystem size is 
increased, implying that the entanglement in the bulk is 
short-ranged. It is this localised nature of the entangle- 
ment entropy around the block edges which gives rise to 



an area law [13, [3 and makes ground states of gapped 
ID systems particuarly amenable to simulation through 
matrix product states (MPS) 

One might expect that boundary effects to the entan- 
glement entropy also have different behaviours for criti- 
cal and non-critical systems, and an interesting question 
then is whether the boundary effects to the entanglement 
entropy are short- or long-ranged, and what the exact be- 
haviour of these are. This is one of the main motivations 
of this work: here we study the entanglement entropy of 
a VBS with general open boundary conditions. We will 
show that, in contrast with the critical XXZ chain, the 
boundary effect to the entanglement entropy in the VBS 
state is short-ranged. We will also show that the satu- 
rated bound for this state is the sum of von Neumann 
entropy of the two boundaries. Furthermore, it is not a 
constant (as is the case for a fixed boundary condition); 
it varies for different boundary conditions. This satu- 
rated bound in ID state corresponds to the area law for 
higher-dimensional systems. 

This model was originally studied by Affleck et al. in 
the context of the Haldane conjecture O) [U) ■ I* ^^^s 
also been the focus of much renewed interest since its gen- 
eralisation to MPS — which have been shown to efhciently 
simulate many ID systems [l^, US] and may be used as 
a variational set in density matrix renormalisation group 
(DMRG) calculations— [2I ill and the discovery that 
its analogue in 2D is a resource for universal quantum 
computation |_23i] . We hope that studying the boundary 
effect on the entanglement entropy will give some further 
insight into this model. Since the boundary effects to 
the block entanglement entropy decays fast, we expect 
the the DMRG method can be applied efficiently to this 
model. 
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II. DEFINITION OF THE VBS STATE 

The spin-1 VBS state with general open boundary con- 
ditions (GOBC) HI takes the form: 



IVBS) = gf 



n 



(44+1-44+1^ 



,k=~Ni + l 



Q?|vac) (1) 



where a\, bj, are bosonic operators, Qf and are respec- 
tively the left and right boundary operators; p,q — ± 
with = a^jv, +1, Qr = ^-jv,+i' Qt = at^w and 



Q~ = b\j^-^ ; I vac) is the vacuum state; and Ni,Nr are 
integer numbers. Since the left and right boundary oper- 
ators are mutually independent, there are altogether four 
different VBS states with GOBC. Note that all sites in 
the spin chain including the left and right boundary sites 
-Ni + 1, L + Nr are spin-l's. Thus this VBS state ^ 
is different from that studied in Ref. Q. We should also 
note one boundary operator, for example, Ql' changes 

the boundary state a^_j^^_f.i\vac) and b^_j^^_f.i\vac) . In 
fact, the state ([T]) is the ground state of the the Hamil- 
tonian studied by Affleck et al. 



H = 



E 

j = -Ni + l 



1 



(s, • s. 



(2) 



III. DIVIDING THE CHAIN 

For convenience in later calculations, we divide this 
ID state into three parts: the left-hand, central and 
right-hand parts. The left-hand part is defined as 
|left,p) = Qf nL-JV,+i(44+i - ^l4+i)|vac). Sim- 
ilarly, the right-hand part is defined as |right, q) = 
I\t=L+i^i"lbl^i - ^la^JQ«|vac). Finally, the central 
part is written |central) = 11^=1(44+1 " 44+i)|vac). 
Note that site 1 appears in both the left and central 
parts, and acts as two spin-l/2's; site L similarly ap- 
pears in both the central and right parts. Thus the whole 
VBS state with GOBC now takes the form | VBS; p, q) ^ 
|left,p)|central)|right,(7). We should note that this is not 
strictly a product state, but that this decomposition is 
valid for our purposes. Double-counting is avoided since 
each bulk spin consists of two spin- 1/2 's and the two 
bosonic operators (spin-l/2's) in one site constitute a 
spin-1 state by Fock space representation. For example, 
terms {a\b\ — 5Qa|)|vac) and {a\b\ — 6|a2)|vac) belong 
to left and central parts, respectively, however, by Fock 
space representation, the product state will create at site 
1 the state (a|)^|vac), (5|)^|vac) and a|6||vac). Thus the 
three parts are connected to constitute the original state 
([1]). Our aim is to now study the von Neumann entropy 
of the reduced density operator of the contiguous spins 
from site 1 to L of the state |VBS;p, q). For this aim, ac- 
cording to the theory of entanglement, the left- and right 



part states |left,p) and | right, g) can be replaced by two 
bipartite states, through the Schmidt decomposition. 

Without loss of generality, we start from the left part 
and consider the entanglement of the quantum state 
|left,p) between site 1 and the rest; i.e. we consider it 
a bipartite state with site 1 as one particle and the rest 
as another particle. According to the Schmidt decompo- 
sition, we can first calculate eigenvalues of the reduced 
density operator of site 1 for state |left,p). 



FIG. 1: Tlie quantum spin chain witli L -\- Ni -\- Nr sites; 
each is spin-1. We calculate the von Neumann entropy of L 
spins in the chain. The whole spin chain is divided into three 
parts: the left, central and right parts. The spin-1 at site 1 
is divided into two spin-l/2's; one in the left part, one in the 
central part. Similarly the spin-1 at site L is split between 
the central- and right parts. 



fc"fe+i - 44+1 ) I vac), we know 



Denote 1*-)^,^+! = (44 
|*-)oi|*-)i2^EL=iI"i)(^ ® CT„i)l^'~)o,2, where 
Wi}i=o is the Pauli group, and we have defined the states 
\ai) — I ® 1^*^). Here (Ti — a\bi + aib\, G2 = 



-ia{bi + iaib{, 0-3 



a\ai 



b\bi and (Tq = a\ai + b\bi, 



the identity. By this result, the state of the left part may 
be written 



|left,p) 



1 



3(7V,-l)/2 



E 



\a-N,+2) 



ao,...,Q!-Nj+2 = l 



-Wi+1,1- 



® \ao) 



(3) 



It is now possible to calculate the site 1 reduced density 
operator. Using the identity J2a=ii^ ®^a)\'^^){^^ 1(1 ® 
aaY — I — |^'^)(5'^| (where / on the l.h.s. and r.h.s. is 
the identity in (D^ and (D^ ® (D^, respectively) we find 



Pi 



= Tri(gf ® /)[-(! - fi)I + m-){^-\]iQf ® ^)^(4) 



where /; = : and the trace is over the first 

Hilbert space. Now we find that the matrix form of 
the reduced density operator of site 1 takes a diag- 
onal form pi = diag(^j'^, ^j~), where we have defined 



= (3 ± /0/3 and fi = (-1/3) 



N,-l 



For different 



boundary operators , state pi is invariant under a 
basis transformation. By entanglement theory, we can 
replace the quantum state of left part by a real bipartite 

state |left,p) \(t)i) = [ J £,^alb{- J £,^bla{ ]\va,c). 
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We find that the reduced density operator pi converges 
to the identity exponentiaUy fast with respect to Ni , and 
thus we can simply consider as a singlet state I^'q i) 
when Ni — > oo. Similarly for right part of the state, we 

have \(t)r) = (^\/C?'^L^L+i ~ |vac), where 

the have a similar definition to the . Thus the 
VBS state in Eq. ([T]) may be rewritten 

L-l 

|VBS) = n (4^1+1 - fcl4+i)l'/'r)|vac), (5) 

k=l 

where indices p, q are suppressed since they do not change 
the result. The validity of the transformation from ([T]) to 
([5|) in studying the von Neumann entropy of contiguous 
L spins can also be checked by a method with matrix 
product state representation introduced in, for example, 
Ref. J2?|. A numerical calculation for small L confirms 
our result. 



LEFT LEFT 




FIG. 2: The quantum state of left part is considered as a bi- 
partite state with one particle in site 1 and the rest as another 
particle. According to the Schmidt decomposition this state 
is equal to a bipartite state \(f>i), and each spin is spin-1/2 at 
site and 1; Finally, the VBS state with GOBC is mapped to 
a state with L + 2 sites with two spin-1/2 boundaries at ends 
and L + 1 



In the case that Ni ^ (x, Nr oo the state takes 
the form |VBS) = nto(4^l+i " ^l4+i)|vac). The en- 
tanglement entropy of this state has been studied pre- 
viously [3]; it was found that there are no boundary ef- 
fects. In this paper, one of our main concerns is to show 
that the VBS state ([1]) does have a boundary effect. Al- 
ready we know that if the block of L contiguous spins in 
Eq. Il]) is far from the two boundary sites the bound- 
ary effect will decay very rapidly (exponentially). We 
now present explicitly the entanglement entropy of these 
L spins. Let us first rewrite the left part state in the 
form 10/) = (V; ® /)|^'~)oi, where the matrix form of VJ 
takes the form V/ = diag(^,'^, ^,~). Note that V) is not 
necessarily unitary. In a similar manner for right part 
state , Vr can also be defined Vr = diag(^~,^+), giving 
= (/ (g) Vr)\'i'^}LX+i- We can then write the VBS 
with GOBC in the form: 

|VBS) = (1^ ® K)o,L+i 51 ® • • • ® 

X((7ai •■•O-at (8)/)|1'")o,L + l, (6) 

where the summation is from 1 to 3 for indices 
ai, ■ • • ,Q!L. According to entanglement theory, the von 
Neumann entropy of the reduced density operator of L 
spins is the same as the von Neumann entropy of the two 
ends (sites and L + 1). The reduced density operator 
of these sites takes the form 



PL^{Vl<» Vr) I -^L + \UVl<S) Vr)^ (7) 



where p = (—1/3)^. 



Expanding this, wc find the following form for the density matrix: 



PL = 




V 





^2^1 










When p is small the eigenvalues of the matrix pi, can be found by a Taylor expansion to be Ai = £,i~S,r (1 ^ p); 
A2 = CrC(l -P); ^3 = CCi^+P) + 0{p^); A4 = CrC(l +P) + Oip^). Recall that \(f>i) is a pure state, so the 
reduced density operators at sites and 1 are the same under unitary transformation (i.e. po = pi); similarly for 
state \4>r)- By checking the eigenvalues of pL, we find that pl = Pq ® Pl+i + 0{p), where the equation is true under 
a unitary transformation. This transformation has no effect on the von Neumann entropy, and thus we suppress it. 

The density operator of pL converges exponentially fast to the tensor product of two ends; the speed is p — (—1/3)^. 

^ 



IV. BLOCK ENTROPY 

Finally we calculate explicitly the von Neumann en- 
tropy of a block of L spins. 



The von Neumann entropy of L spins converges to the 
von Neumann entropy of two ends exponentially fast with 
convergence p = (—1/3)^. The exact form of 0{p) is 
written below [31.]. For a pure bipartite state \4'i), we 



S{pL)^S{po) + S{pL+i) + 0{p). (8) 



4 



3-fi 
6 



3±Aiog 



3+/1 

6 ■ 



f 

and similarly 



know S{po) = S{pi) = -^log 

When /; is small, we find S{po) ^ 1 
for S{pl+i)- So when the distances between the block of 
L contiguous spins and the two ends are large, we know 
that 



S{p~l) 
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+ 0{p) + 0{fr) + 0{fr). (9) 



Note that when p is comparable with fi and /,., then 
0{p) may have contributions for both /; and fr- Consid- 
ering all of these points, we conclude that boundary ef- 
fects decay exponentially when the distances between the 
subsystem and the boundaries increase. In the case that 
there are no boundary operators, the trace in Eq. ^ is 
over the identity, and pi reduces to /. Hence the bound- 
ary effects to the entanglement entropy never arise; this 
is the case previously studied 0] ■ We remark for ID VBS 
with GOBC, the terms corresponding to the topological 
entropy in Ref. [2^, [2^ are the boundary terms appeared 
in Eq.®, and they remain unchanged for L —>■ 00. 



V. COMPARISON WITH A CRITICAL SYSTEM 

It was found numerically that for spin-1/2 XXZ 
chains the entropy takes the form S{L, N) — Su{L, N) + 
{—1)^Sa{L,N), where the second term arises from the 
boundary conditions, and can be written SAiL,N) — 
l/( sm{2nNn.r/N)N/TT)'^ . Note that the notation of Ref. 
[15| has been changed slightly here and N = Nr + Ni+ L, 
Ni = Nr = Nn^r] K depends on the anisotropy parameter 
in the XXZ chains and i^T = 1 for an XX chain. We 
find that the boundary term decays slowly and is quasi- 
long-ranged, while for the VBS studied in this paper, the 
boundary terms in Eq. © decay exponentially fast. 



VI. TWO-SITE ENTANGLEMENT BY 
NEGATIVITY AND REALIGNMENT 
CALCULATIONS 

We next study the boundary effects to entanglement 
between only two spins in the bulk. Consider spins 1, L. 
The previous method still works; i.e. we can reduce the 
length of the chain from Ni + L + Nr sites in Eq. ([T]) to 
a chain with only L + 2 spins in Eq. ^ . This transfor- 
mation does not change the entanglement between two 
spins at sites 1 and L. We find that the density operator 
of the bipartite state is written 

Pi,L =^|aiaL)(aiaL|Tr(V^0V;)(cr„i cr^J 



l-p 



I + p\'^ ){^' 



Ki )t(i^®i/,)t (10) 



The explicit form of this density matrix is complicated. 
However since the boundary matrices Vi and Vr become 
identities exponentially with Ni and Nr, we know that 





Nr = 1 


TV,. = 2 


Nr=3 


Nr ^4 


Nr = 00 


Ni = 1 


1.45919 


1.50111 


1.43456 


1.45142 


1.44170 


Ni=2 


1.50111 


1.05433 


1.15504 


1.11552 


1.12486 


Ni=3. 


1.43456 


1.15504 


1.00609 


1.05018 


1.03861 


Ni=4 


1.45142 


1.11552 


1.05018 


1.00068 


1.01252 


N, = oo 


1.44670 


1.12486 


1.03861 


1.01252 


1 (exact) 



TABLE I: Nearest-neighbour negativity 





Nr = 1 


iV,. = 2 


Nr=3 


Nr ^4 


iV,. = 00 


N = 1 


0.37393 


0.23445 


0.19692 


0.20032 


0.19861 


Ni^2 


0.23445 


0.03974 


0.02631 


0.02194 


0.02221 


Ni=3 


0.19692 


0.02631 


0.00439 


0.00293 


0.00247 


Ni^4 


0.20032 


0.02194 


0.00293 


0.00049 


0.00027 


N = 00 


0.19861 


0.02221 


0.00247 


0.00027 


(exact) 



TABLE II: Nearest-neighbour entanglement by realigment 



the boundary effects to this density operator decay ex- 
ponentially if the separation of these two spins with the 
boundaries increase. Thus the case that iV^ or Ni is small 
already can provide enough information about the en- 
tanglement between the considered two spins. The cases 
Nr 00 and iV; — s- 00 are special since they correspond 
to the cases that there are no boundary operators 
and , respectively. 

In order to quantify the two-spin entanglement in this 
case, we shall use negativity, M |23|. We find that there is 
no entanglement for non-nearest neighbouring spins (re- 
gardless of boundary) . The nearest-neighbour negativity 
for a range of TV; , Nr is presented in Table HI where a 
factor 1/9 is omitted. 

Curiously, we see that while boundary opera- 
tors decrease the block von Neumann entropy — see 
Eq.® — they increase the nearest- neighbour negativity. 
Roughly, this can be understood that the entanglement 
is monogamous [1^ [s^l, i.e., it can not be shared freely 
by many parties. A simple example about the monogamy 
of entanglement is that, suppose A, B and C are three 
parties, if A and B are maximally entangled, A and C 
will be separable. The result in this paper shows that 
the boundary operators have effect on the entanglement 
sharing in this many-body system. 

In the table above, we see that for finite Ni^r, it is 
always that case that M > 1/9. Asymmetric boundary 
operators may also increase the entanglement. For exam- 
ple, M = 1.45919/9 for Ni ^ Ni ^ 1, but TV = 1.50111/9 
for Ni = 1, Nr — 2. Now Af does not vary monotonically 
with Ni^r', e-g- when Nr varies from 1 to 00 (for Ni = 1), 
J\f oscillates. This seemingly new result may potentially 
be useful: nearest-neighbour entanglement may be con- 
trolled by tweaking boundaries. 

The drawback of the negativity is that the bound en- 
tanglement cannot be detected and quantified. A com- 
plementary quantity derived from the realignment sep- 
arability criterion can partially solve this problem [28| . 
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For a two-site density operator p, this quantity is defined 
as = (||i?(/r?)| | — 1)/2, where the matrix R{p) is obtained 
from the density operator p by the reaUgnment method, 
and II • II is the trace norm. The larger of the quanti- 
ties M and TZ gives a lower bound of the concurrence C 
for a mixed state in arbitrary dimensional systems; i.e. 
C > max{A/',7e} 

Using this measure, we still find zero entanglement 
between detected for non-nearest neighbouring spins re- 
gardless of boundary. The entanglement of nearest neigh- 
bouring spins by realignment is presented in Table HIl and 
a factor 1/9 is also omitted. We also still find that the 
boundary operators increase entanglement TZ, and indeed 
no entanglement is found for the case without boundary 
operators: this is due to the limitation of the realign- 
ment method. Finally let us remark that for the model 
studied in this paper, the negativity provides a stronger 
lower bound for the concurrence. Since the concurrence 
for a general mixed state in (g) is difhcult to find, 
the entanglement measures by negativity and realignment 
are widely accepted. 

VII. CONCLUSION 

In summary, we studied the boundary effects to the en- 
tanglement entropy and the two-site entanglement for the 



spin-1 VBS state which is the ground state for a gapped 
model. We showed that the boundary effects are short- 
ranged, i.e. they decays exponentially in the distance be- 
tween the subsystem considered and the boundary sites. 
This is different from the case of XXZ chain which is 
a gapless model. The two-site entanglement was stud- 
ied by two entanglement measures, the negativity and 
the realignment method. For the VBS state, we find the 
boundary operators decrease the block von Neumman 
entropy but increase the the nearest-neighbour entangle- 
ment measured by negativity and realignment. 
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